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ABSTRACT 

NASA has been studying future transport concepts, envi- 
sioned to be technically realizable in the timeframe of 2020- 
2030, to meet environmental and performance goals. One con- 
cept receiving considerable interest involves a propulsion system 
embedded into a hybrid wingbody aircraft. While offering signif- 
icant advantages in fuel savings and noise reduction by this con- 
cept, there are several technical challenges that are not encoun- 
tered in the current fleet and must be overcome so as to deliver 
target performance and operability. One of these challenges is 
associated with an inlet system that ingests a significantly thick 
boundary layer, developing along the wingbody surface, into a 
serpentine diffuser before the flow meeting fan blades. The flow 
is subject to considerable total pressure loss and distorted at the 
fan face, much more significantly than in the inlet system of con- 
ventional aircraft. In our previous studies [12], we have shown 
that through innovative design changes on the airframe surface, 
it is possible to simultaneously increase total pressure recovery 
and decrease distortion in the flow, without resorting to conven- 
tional penalty -ridden flow control concepts, such as vortex gen- 
erator or boundary layer bleeding/suction. In the current study, 
we are interested in understanding the following issues: how the 
embedded propulsion system performs under a crosswind condi- 
tion by studying in detail the flow characteristics of two inlets, 
the baseline and another optimized previously under the cruise 
condition. With the insight, it is hoped that it can help in the 
follow-on study by devising effective strategies to minimize flow 
distortion arising from the integration of an embedded- engine 
system into an airframe to the level acceptable to the operation 


of engine fan. 

NOMENCLATURE 

a speed of sound. 

A Area. 

D Diameter. 

DPCP Circumferential distortion descriptor recommended by 
SAE. 

L Length measured in the x-direction. 

M Mach number. 

P t Total pressure. 

Re Reynolds number, 
p Density. 

x,y,z Axial (horizontal), spanwise and vertical direction. 
Subscripts 

0, °o Far upstream or “free” stream. 

2 Aerodynamic Interface Plane (AIP) of engine. 
c Inlet lip. 
i Throat. 


INTRODUCTION 

NASA, under the current Fundamental Aeronautics Program 
(FAP), is considering several commercial transport configura- 
tions beyond current aircraft in the Subsonic Fixed Wing (SFW) 
Project m , with objectives to achieve considerable improvement 
in performance and environmental impacts. Table [I] summarizes 
the specific goal metrics targeted to reduce noise, NOx emission, 
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and fuel consumption before 2030. To achieve these demanding 
goals, non-conventional concepts are called for; but technology 
gap is too big that it requires evolutionary approach by focusing 
various concepts and technologies needed in the next three gen- 
erations of aircraft, respectively named as N+l, N+2, and N+3. 
Noticeably, considerable reduction in each category of |T] is re- 
quired in N+2 (relative to Boeing 777-200 and GE90 engines) 
and N+3 (relative to Boeing 737-800 and CFM56-7B engines). 
In this study, concepts for N+2 is our interest. A concept that has 
potential to achieve these metrics and has been under intensive 
study is the hybrid wing body (HWB) airframe with a tightly in- 
tegrated propulsion system, see Fig. Q] The inlet is non-circular 
at the entrance and the entering flow, no longer uniform or free 
of disturbances, and is now carrying with it a boundary layer de- 
veloping along the fuselage; the inlet is thus known as boundary- 
layer-ingesting (BFI) inlet. 


TABLE 1 : NASA’S SUBSONIC FIXED WING PROJECT. 


TECHNOLOGY 

TECHNOLOGY GENERATIONS 
(Technology Readiness Level = 4-6) 

BENEFITS* 

N+1 (2015) 

N+2 (2020**) 

N+3 (2025) 

Noise 

(cum margin rel. to Stage 4) 

-32 dB 

-42 dB 

-71 dB 

LTO NOx Emissions 
(rel. to CAEP 6) 

-60% 

-75% 

-80% 

Cruise NOx Emissions 
(rel. to 2005 best in class) 

-55% 

-70% 

-80% 

Aircraft Fuel/Energy Consumption* 
(rel. to 2005 best in class) 

-33% 

-50% 

-60% 


* Projected benefits once technologies are matured and implemented by industry. Benefits vary by vehicle size and mission. N+1 and N+3 values 
are referenced to a 737-800 with CFM56-7B engines, N+2 values are referenced to a 777-200 with GE90 engines 
** ERA'S time-phased approach includes advancing "long-pole" technologies to TRL 6 by 2015 
t CO 2 emission benefits dependent on life-cycle C0 2e per MJ for fuel and/or energy source used 

The concept of ingesting boundary layer into a propulsion 
system is not entirely new, originated in 1947 by Smith and 
Robert (U, who concluded the performance advantages with a 
significant longer range by submerging turbojet engines into the 
wing or fuselage. Subsequently further analyses were carried 
out at the Douglas Aircraft Company to assess the benefits of 
boundary layer ingestion to propulsion efficiency OH . In 2002, 
Eiebeck m chronicled the development of the blended wing 
body, in which the boundary-layer-ingestion propulsion system, 
assuming it was possible to provide a uniform flow and efficient 
pressure recovery at the fan face. Because of recent demands 
in reducing noise and fuel consumptions, renewed interest in 
embedded propulsion systems has arisen, significant efforts and 
investigations have been invested by various institutions, such 
as NASA’s SFW project, Cambridge-MIT’s Silent Aircraft Ini- 
tiative D, and European Commission’s New Aircraft Concepts 
Research (NACRE) 0 and many reports have been produced, 
e.g., mmwrnmm, documenting benefits and technology 
barriers and risks before the concept is realized in future com- 


mercial transport. Especially, rich depository of studies reports 
can be found from the NASA and SAI, in which various specifi- 
cations are put forth and configurations published. 

To fully reap the advantages of a B LI propulsion system, 
it must be placed close to the airframe so that no discernible 
wakes arise individually from them. Thus, it requires a geometri- 
cal change from being non-circular transitioning to circular and 
offsetting of its centerline, as shown in Fig. Q] and O which is 
designated as N2B concept. As depicted in Fig. [lj the configura- 
tion is designed to have three separate engine systems, with the 
center one situated further downstream from the side ones. Each 
nacelle will house a turbojet engine in the center flowpath, while 
the two side flowpaths will only include respectively a fan driven 
by the core engine, as suggested in Fig. |2l 

As a result of mutual interference between flows over the air- 
frame and nacelle and entering the inlet, the flowfield becomes 
enormously complex; asides from being turbulent, it typically 
involves flow separation, 3D vortical flow and possibly unsteadi- 
ness, in the inlet duct. This results in an increased nonuniformity, 
total pressure loss and distortion, which lead to additional con- 
cerns about blade mechanical loads, fan stall and even engine 
surge, and engine noise. Moreover, because the design of the 
BLI system is targeted at cruise conditions, new aerodynamic 
characteristics at off-design conditions arise and may have more 
significant impacts than the conventional aircraft configuration. 
A particular concern is about the aerodynamic performance in 
the BLI inlet when under the influence of crosswind at low speed 
conditions during takeoff and landing. This is the subject of the 
present study. 



FIGURE 1 : N2B CONCEPT CONFIGURATION (3). 


Past studies of flow in BLI inlet include experimental and 
computational works, such as fl5l and |d6, 17 ] for inline and 
offset (S-shape) inlets, in which various flow control concepts 
were explored so as to reduce flow distortion. With consider- 
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FIGURE 2: PROPULSION SYSTEMS AND NACELLE CON- 
FIGURATION IN N2B. 


able amount of ingested boundary layer (typically 30% of in- 
let height), distortion at aerodynamic interface plane (AIP) in- 
evitably exceeds the level tolerable by the engine fan, thereby ne- 
cessitating controlling, even reducing of boundary layer growth. 
Various flow control concepts have been studied for inlet applica- 
tions; most chose vortex generators (VGs) [ 18, 15] 19 ], and some 
used active control 120111 11 or bleeding. However, they all incur 
losses and costs, such as depriving of useful air supply, requir- 
ing additional power, keeping intact of VGs and bleed holes, etc. 
Moreover, determining the location, number, size and shape of 
VGs to be employed for achieving optimal results and maintain- 
ing the same level of effectiveness for a wide range of operat- 
ing conditions has proven to be difficult because the number of 
design variables is large and the flow is extremely complicated. 
A systematic approach to quantitatively defining the above fac- 
tors is perhaps the only robust route. The design of experiment 
method has been used in OiiQj ] . Recently, genetic algorithms 
(GAs) together with data mining were employed to determine 
the locations and orientations of vortex generators so as to effec- 
tively reduce distortion and increase total pressure recovery ED- 
While GAs have several advantageous features, including ease 
of implementation and allowing searching for a global optimal 
point, they are computationally expensive for a large problem. 
In another study m, we chose the adjoint sensitivity approach 
for minimizing distortion by redesigning the airframe surface to 
“condition” the approaching boundary layer to the inlet. With 
this new intriguing surface shape, not only is the distortion at the 
AIP significantly reduced, but also the total pressure recovery in- 
creased considerably. These benefits are realized only through 
a fundamental change in flow structures, without incurring addi- 


tional mechanical devices or loss of useable compressed air for 
boundary layer bleeding. 

With the operability of using the BLI concept demonstrated 
for the design condition, we now shall focus on understanding the 
flow characteristics at off-design conditions, of special impor- 
tance is to assess the performance under crosswind at low speed 
flight during takeoff and landing. This subject has received much 
less attention than that for the cruise condition, rightly so because 
an aircraft spends its vast majority of time in cruising. However 
it is hugely important to ensure safety during this very brief, yet 
critical, phase of the flight. Safe landing and takeoff may be re- 
quired with crosswind speed on the order of 15-20 m/s (roughy 
30-40 knots) and the associated sideslip angle may be of 15 to 20 
degrees for conventional transport. Recently, several detailed ex- 
perimental and computational analyses have been reported, see 
e.g. E2I23E31; they all focus on an isolated nacelle at ’’zero” 
cruise speed and subject to crosswind. Severe flow separations 
inside and outside the nacelle, resulting in significant distortion 
to the flow before AIP, are observed. Considerable hysteresis is 
encountered as the crosswind changes angles or magnitude l23l . 
The deterioration in performance is expected to be even worse in 
the case of BLI inlet. It is unclear how the performance deteri- 
orates with respect to relevant factors, such as crosswind angle, 
magnitude, flight speed, etc. Thus, it is an even more important 
step to make a definitive quantitative analysis and to devise feasi- 
ble practical changes to the embedded nacelles/inlet and the air- 
frame surface surrounding it. In this paper, we shall take the first 
step to characterize the resulting flow phenomena and assess its 
aerodynamic implications at the AIP so that these findings will 
lay the foundation and provide the direction for the subsequent 
optimal redesign. 

The remainder of the present paper will cover the following 
topics: (1) Choice of CFD requirements for accurately evaluat- 
ing the complex flowfield in the BLI system, (2) Statement of the 
problem considered in this study (3) Validation of the CFD ca- 
pability for capturing key aerodynamic performance metric, (4) 
Aerodynamic performances of the original and optimal designs 
at the cruise condition, (5) Characterization of aerodynamics un- 
der off-design crosswind condition, and (5) Future research. 

CHOICE OF CFD REQUIREMENTS 

The N2B configuration is designed to fly at high subsonic 
speed and the boundary layer ingested into the inlet is turbu- 
lent. The CFD capability needed for reliably handling this type 
of flow should possess two main technical elements: (1) reliabil- 
ity for accurately predicting complex turbulent flow from low to 
cruise speeds, and (2) flexibility for resolving different geome- 
try topologies and viscous layers without sacrificing grid quality 
(orthogonality and mesh spacing). The analysis is based on solv- 
ing three-dimensional Reynolds-averaged Navier-Stokes equa- 
tion with closure by the two-equation K — co - SST model ESI 
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The basic CFD code has been validated for many problems in 
the past, as documented in the PhD dissertation by Lee [26]; it 
has been further augmented and comparison with experimental 
data for the boundary layer ingestion inlet has been verified in our 
previous study m. For the present study with emphasis on cross- 
wind condition, low Mach flow dominates in regions internal and 
external to the nacelle, for which the conventional CFD schemes 
are afflicted with significant numerical dissipation, thereby con- 
tributing to inaccuracy and deterioration in convergence. An ad- 
vanced upwind method, AUSM + -up l27l . capable of handling 
flow from low to high speed regimes, has been implemented in 
the code. The time discretization is accomplished by the implicit 
LU-SGS scheme [28]. Since the k — co-SST model is regarded 
as the state of the art model for a wide range of flow phenom- 
ena, no attempt in this study to consider other models or other 
formulations such as large-eddy simulation, as this subject mat- 
ter is beyond the scope of the current investigation. It is noted 
that the variables presented are non-dimensionalized, unless in- 
dicated otherwise; the density, velocity and pressure are non- 
dimensionalized by the freestream quantities, Poo,<2oo, and po^cioo. 

The need for flexibly handing complex geometry can be re- 
alized by utilizing the overset grid strategy m, in which two or 
more grid blocks overlap and each has an independently gener- 
ated grid, thereby allowing a better control of grid quality and 
flexibility for adapting to geometry. 


PROBLEM STATEMENT: CFD SOLUTION OF BLI INLET 
FLOW 



FIGURE 3: SCHEMATIC OF THE BASELINE BLI INLET 
UNDER STUDY: A SIDE VIEW SHOWING INCOMING 
BOUNDARY LAYER AND PRIMARY GEOMETRICAL NO- 
TATION. ALL DIMENSIONS ARE IN INCHES, DRAWINGS 
NOT TO SCALE, X-COORDINATE IS ALIGNED WITH THE 
SYMMETRY PLANE AND X=0 COINCIDES WITH THE IN- 
LET HIGHLIGHT. C3 


nated as type A in mn this inlet has been extensively studied 
experimentally and numerically. Its geometry is schematically 
shown in two views in Fig. [3] together with relevant dimensions, 
where stations c, /, and 2 respectively referring to the lip (high- 
light), throat, and AIP of the inlet. The X-Z plane coincides with 
the inlet symmetry plane with the X-coordinate pointing down- 
stream. the inlet highlight is defined at X = 0. The key geometric 
parameters are: the diameter at AIP D 2 = 2.448 in, the throat is 
located at X[ / D 2 = 0.1949 and the AIP at X 2 /D 2 = 3.339, the 
throat height Hi/D 2 = 0.6957, and the area ratio A 2 /A; = 1.070. 
The offset of the inlet is 0.77357)2. The figure also illustrates 
schematically an incoming boundary layer to the inlet, with an 
approximate thickness of 30% of entrance height. We note that 
in most line plots we interchangeably use D for D 2 . 



FIGURE 4: COMPUTATION DOMAIN USED IN THIS 

STUDY AND AN ENLARGED VIEW OF THE INLET 
ENTRANCE, SHOWING THE OVERSET GRID SYSTEM, 
EACH COLOR DENOTING AN INDIVIDUAL BLOCK 
GRID. 


To encompass flow domains inside and outside the inlet and 
allow accurate resolution of viscous layers, using overset grids 
has advantages in generating quality grid and dealing with dif- 
ferent topologies with ease. The former becomes especially at- 
tractive in shape optimization since the grid will be changing 
throughout the optimization process. Figure [4] shows the com- 
putation domain used in this study and an enlarged view near 
the BLI inlet entrance, depicting the overset mesh system, con- 


In this study, we consider a circular type BLI inlet (desig- 


1 While the inlet entrance in N2B is more of super-elliptic shape, it is believed 
that qualitative phenomena in both cases are similar and more importantly the 
N2B inlet is still a subject of conceptual study and there is no data available yet 
for validation. 
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taining 16 grid blocks, 7 configuration blocks and 9 background 
blocks, in which individual grid blocks are designed to resolve 
the curvature of inlet lip, the “inviscid” region outside the inlet, 
the core region inside the inlet background, and the boundary 
layer near the inlet entrance, the viscous region next to the inlet 
wall. The entire mesh has about 10.7 millions grid points in total 
and more than doubles the grid used in the previous study where 
only a half domain is needed by employing symmetry since no 
crosswind is involved. 

As the main interest is to investigate how well the BLI inlet 
performs under the crosswind condition. The condition consid- 
ered critical is during takeoff and landing phase. The crosswind 
strength is specified to be 40 knots (20.78 m/s) and the sideslip 
angle is 30 degrees. At the sea level dry air condition, this yields 
a crosswind Mach number of 0.06 and freestream Mach number 
of 0. 1212. The boundary-layer thickness at the entrance of the in- 
let is targeted to be 0.35 H c , an expected value in a realistic flight 
situation. ® Short of knowing the boundary-layer flow on actual 
configuration, we assumed that a turbulent flow is developing on 
a flat plate at a distance of 20£>2 upstream of the inlet, such that 
the required boundary-layer thickness is obtained. Since the flow 
at AIP is subsonic, as required for the fan operation, the usual 
strategy of imposing static pressure there is used. However, mass 
flow rate is the parameter customarily known in engine opera- 
tion, instead of static pressure. This can be bridged by assuming 
a locally ID isentropic condition, yielding a direct relationship 
between mass flow rate and static pressure (e.g., see II). 

The outflow condition is specified at a somewhat arbitrary 
interface plane with the fan, if the interaction between inlet and 
fan is neglected in the analysis. The computational condition at 
this “virtual” AIP can be handled several ways, such as specify- 
ing a mass flow rate or static pressure. Since there is no known 
information about how the static pressure distribution would be 
radially, it is important to minimize influence of any a priori 
specification of its value and profile. Hence, the computation 
domain is extended beyond AIP with additional straight section 
of same diameter and of length La = (X 2 — 2Q). It is verified af- 
ter the solution is converged that in this extended section the flow 
becomes nearly unchanged towards the end and the static pres- 
sure becomes uniform, but the pressure at AIP still has radial as 
well as circumferential variations. 

VALIDATION AND COMPARISON WITH DATA AT DE- 
SIGN CONDITION 

To ensure the predictive capability for dealing with complex- 
ities relevant to the problem of our interest, we have carried out 
extensive validation tests and the results were documented previ- 


2 It is assumed that there is no spanwise variation in the flow approaching the 
inlet since the inlet is attached to an infinite plate upstream. In actual aircraft, the 
wingbody will have a spanwise variation in shape, thus adding another dimension 
of complexity to the flow. 


TABLE 2: COMPARISON OF COMPUTED INLET PERFOR- 
MANCE WITH MEASURED DATA. 


Reynolds number 

Total pressure 
recovery ratio 

Distortion 

Mass flow ratio 

A 0 / A c 

3.8xl0 6 O 

0.952 

0.0540 

0.534 

3.8xl0 6 present CFD 

0.956 

0.0596 

0.533 

2.2x1 0 6 present CFD 

0.943 

0.0630 

0.527 


ously in Q. From Table |2l it is evident that the employed CFD 
code is capable of providing a reliable solution insofar as the in- 
let performance metric is concerned. A lower Reynolds number 
condition was also computed because it was suggested in a pre- 
vious study by Berrier and Allan l30ll that this lower value gave 
a closer agreement in pressure distribution on the bottom wall. 
However, this table indicates that (1) lowering Reynolds number 
produces less total pressure recovery, higher distortion and lower 
mass flow ratio, this higher-loss phenomenon suggests an exis- 
tence of a larger separated region, as manifested by the fact of a 
flatter pressure rise. However, with this choice of lower Reynolds 
number, the inlet metric deviates noticeably from the data. Hence 
we conclude that the discrepancy in surface pressure distribution 
might have to be attributed to other factors. 

OPTIMIZATION ON CRUISE CONDITION 

Shape optimization of the inlet bottom wall has been car- 
ried out in our previous studies (HE) based on the adjoint sensi- 
tivity formulation to minimize flow distortion DPCP avg at AIP 
on cruise condition: Moo = 0.85 , Reoo = 3.8vl0 6 . DPCP avg is 
a circumferential distortion descriptor, defined according to the 
SAE standard ED by averaging total pressure distortions over 5 
equally area- weighted rings: 



where P t j represents the average over the 12 equally spaced 
points on ring i and Pt,i low is the average of total pressures be- 
low P t i. 

Here we include representative and revealing results, which 
demonstrate a considerable simultaneous improvement in distor- 
tion reduction and total pressure recovery. The interested reader 
is encouraged to find details in the cited references. 

Figure[5]displays the “oil-flow” patterns of both designs; the 
original design gives rise to a significant “push” of streamlines 
from sidewalls, because of low-momentum-fluid blockage, while 
the streamlines in the optimal design follow the contours of in- 
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let. These are a result of making a series of “sinusoidal-like” of 
varying amplitudes and wave lengths modulations on the bottom 
wall, in a patch defined by 

-1.8 < — <0.5 and — 1.0 < — <1.0 (2) 

~ D 2 ~ ~ D w ~ 

where D w (x) is the width of crosssection at location x 


Baseline 


Optimized 



FIGURE 5: “OIL FLOW” PATTERN OF THE ORIGINAL 
(TOP) AND OPTIMAL (BOTTOM) INLETS, SUPERIM- 
POSED ON PRESSURE CONTOURS. EH 


Figure [6] shows that the optimized design has consistently in 
all off-design conditions far smaller distortion values, 50 to 70% 
less than the baseline values. Also notice that at the pressure ratio 
Pb/Pt, o = 0.9353 the baseline model approaches unstart with a 
much reduced mass flow ratio, but the optimized model appears 
still functioning. Furthermore, the total pressure contours at AIP 
shown for two mass flow ratios also confirm the superiority in 
pressure recovery after optimization. 



Ao /A, 

FIGURE 6: COMPARISON OF FLOW DISTORTIONS AT 
OFF-DESIGN CONDITIONS. EH 


In what follows, we present the resulting performance of the 
inlet in terms of DPCP avg , total pressure recovery (P t /P t , «,), and 
flow rate, but at takeoff low speed with crosswind. 


CROSSWIND CONDITION: CHARACTERIZATION OF 
AERODYNAMIC PERFORMANCE 

As seen above, redesigning the airframe surface near and 
ahead of the inlet entrance has achieved the goal of improv- 
ing the inlet aerodynamic performance to the extent tolerable to 
the fan. It is now to ensure that the BLI inlet deliver needed 
performance under off-design conditions, especially with regard 
to safety under severe conditions, such as takeoff and landing 
with crosswind. As indicated earlier, the condition considered 
is: cross wind Mach number of 0.1 and sideslip angle of 30 de- 
grees, resulting in Reoo ,£> 2 = 0.58vl0 6 . For this analysis, we can 
no longer use a half plane by assuming symmetry, the complete 
nacelle/inlet must be considered, see the overset grid system in 
Fig. [4] The back pressure ratio is set as Pb/P OQ = 0.7373. The ap- 
proaching boundary layer is generated to give a thickness of 33% 
of highlight height at the inlet entrance. Fine grid spacings are 
ensured near a solid surface in each grid block so that viscous- 
dominated regions are adequately resolved; a representative wall 
distance of a first grid point measured by y + is on the order of 
unity, specifically at the inlet entrance it is less than one. 

In what follows, we shall not only focus on the performance 
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of the baseline BLI inlet and the underlying flow characteristics, 
but will also accompany with those from the inlet optimized for 
the cruise condition flUEl, so that we shall give a quantitative 
assessment as to whether the previous optimization holds its ad- 
vantages over the baseline geometry under the current crosswind 
condition. For the sake of conveniently referring to this geom- 
etry, we shall here afterwards use the term “optimized” inlet to 
mean the geometry previously optimized for the cruise condi- 
tion mm, but it should be emphasized that it is not optimized 
for the current crosswind condition. 



Free Stream Direction 


FIGURE 7: LOCATIONS OF ’’NUMERICAL” PRESSURE 
RAKES FOR SHOWING PRESSURE DISTRIBUTIONS IN 
THE FOLLOWING FIGURE. 


First, we present the surface pressure distributions on four 
representative circumferential sections, as designated in Fig. [7 ] 
on both internal and external surfaces of the inlet. The surface 
between lines 2 and 3 will be denoted as leeward side and the 
opposite (between lines 1 and 3) is the windward side; the sur- 
face 1-4-2 is the bottom surface of the inlet, which eventually 
evolves into the lower half of the circular duct. It is evident that 
there are considerable variations in the circumferential direction, 
between the internal and external surfaces. The pressure distribu- 
tions on the internal and external surfaces on these four sections 
are displayed in Figs. [8] and [9] respectively. On the windward 
side (section 1) the external surface has a stagnation point down- 
stream the nacelle highlight, the pressure there, see top figure 
in Fig. 0 is only slightly greater than EL because the dynamic 
pressure is small due to a low speed, and it reduces to a constant 
value of essentially EL since the outer surface has no changes 
in slope or curvature. While entering the internal inlet, the flow 
however undergoes a rapid expansion as it negotiates around the 
lip and continues to “diffuse” downstream, accompanying with 



x/D 


Inlet 



(x/D-3.3) 



t 



FIGURE 8: PRESSURE DISTRIBUTIONS ON CIRCUM- 
FERENTIAL SECTIONS, ALONG LINES 1 AND 2 (OP- 
POSING WINDWARD AND LEEWARD INLET/AIRFRAME 
JUNCTURES), ON INTERNAL AND EXTERNAL (HAVING 
HIGHER PRESSURE VALUES) INLET SURFACES. 


increasing pressure, which could lead to flow separation (as will 
become evident later in Fig. ns. The external surface on the lee- 
ward side (section 2) yields a similar behavior as in the windward 
side, see bottom figure in Fig. 0 but the internal surface sees a 
much reduced expansion in pressure and the pressure rise down- 
stream is milder, thereby less likely giving rise to flow separation. 
Also included in both plots are the results from the optimized in- 
let. The trends of pressure distributions between the baseline and 
optimized inlets in these two sections are very close, except a 
somewhat reduced pressure expansion in the windward side of 
the optimized inlet. 
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Figure|9]reveals that the pressure distributions on the internal 
surfaces of the top section (section 3) and bottom section (section 
4) are more complex than those just discussed; on section 3 we 
see an extended pressure decrease following the initial rapid rise 
near the throat and a final rise to the specified AIP pressure, un- 
like those observed in sections 1 and 2. Again, the results of 
both baseline and optimized inlets are nearly indistinguishable, 
as expected since the modification in geometry is made on the 
bottom wall only. However, the geometry effect is now observed 
on the bottom wall pressure (section 4), the effect remains well 
into the downstream region of the inlet, up to X/D 2 ~ 0.5, where 


inlet 



J _l 1 I I L 

0 12 3 4 

x/D 



FIGURE 9: PRESSURE DISTRIBUTIONS ON CIRCUMFER- 
ENTIAL SECTIONS, ALONG LINES 3 AND 4 (OPPOSING 
TOP AND BOTTOM WALLS), ON INTERNAL AND EXTER- 
NAL (HAVING HIGHER PRESSURE VALUES) INLET SUR- 
FACES. 





p: 0.30 0.35 0.40 0.45 0.50 0.54 0.59 0.64 0.69 

FIGURE 10: COMPARISON OF PRESSURE CONTOURS ON 
THE SYMMETRY PLANE INSIDE THE BASELINE AND 
OPTIMIZED INLETS. 


the modification of wall shape terminates. The flow here behaves 
differently from that on section 3; the pressure has a milder ex- 
pansion initially and then steadily climbs till the AIP. 

Next in Figs. [lOl we show the pressure contours on the sym- 
metry plane; the flow is accelerating primarily inside the inlet for 
this flow configuration, as also revealed in previous figure Fig. [9] 
and the Mach number contours (not shown) remains by and large 
above 0.6, only around the second bend on the bottom wall (near 
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AIP) is its value lowered. The contours between the baseline and 
optimized inlets are almost the same, except on the bottom wall 
near the entrance where the wall is reshaped. 


Baseline Optimized 



0.850 0.874 0.897 0.921 0.945 0.968 0.992 


FIGURE 11: COMPARISON OF PRESSURE CONTOURS 
ATAIP OF THE BASELINE AND OPTIMIZED INLETS. 


Figure Qj] compares the total pressure contours at AIP of 
baseline and optimized inlets; both plots clearly exhibit a great 
deal of asymmetry, thereby contributing distortion, and have a 
low pressure region skewed towards the windward side of the 
nacelle. The pressure level in the optimized inlet however has 
lower value, suggesting a lower pressure recovery and higher dis- 
tortion. The origin of this low pressure region can be traced to 
the flow entering the inlet, as shown in Fig. |T2j where Mach 
number contours are displayed at several streamwise sections 
superimposed with particle traces inside the inlet. A pocket of 
low momentum fluid was observed in the optimized inlet, which 
in turn contributes to the thickening and diffusion of low mo- 
mentum boundary layer fluid, thus losing total pressure along 
the way. 

Table [3] summarized the performance metrics of both inlets, 
it shows that the optimized shape, as is, does not hold advantages 
under crosswind condition. However, this does not necessarily 
mean that optimization fails under this condition, it may require 
considering an extended design space, in our case it is the area 
upstream of inlet entrance in the spanwise direction. Neverthe- 
less, the pressure recovery is only slightly deteriorated. 

Examination of the velocity field at various cross-sections 
reveals interesting characteristics of the flow, the results of both 


TABLE 3: COMPARISON OF PERFORMANCE BETWEEN 
BASELINE AND OPTIMIZED INLETS. 


Inlet 

Total pressure 

Distortion 


recovery ratio 


Baseline 

0.963 

0.0169 

Optimized 

0.958 

0.0515 


baseline and optimized inlet are represented in Fig. \m The 
flow is clearly asymmetrical with respect to the symmetry plane, 
due to the crosswind component. The cross-sectional velocity 
components (projection of the velocity vector onto a x=constant 
plane, leaving only the y and z components) clearly portraits a 
flow from the windward to leeward sides in the exterior of the 
nacelle; it also exhibits that a separated flow region, denoted as 
the front separation, first appears only in the windward stagnation 
juncture, see the plot at X/D 2 = 0.5. Meanwhile a second sep- 
arated region is beginning to form, visible at X/D 2 = 2.5 on the 
top of the nacelle, roughly 135 degrees from the front stagnation 
point, and is growing in size, as clearly displayed at X/D 2 = 3.3. 
This second separated region does not appear connected to the 
front separation in the sense that it is not morphed from the lat- 
ter. Also its size eventually decrease to zero further downstream 



Mach: 0.05 0.16 0.27 0.38 0.49 0.60 0.71 0.82 

FIGURE 12: COMPARISON OF CROSS-SECTIONAL MACH 
NUMBER CONTOURS BETWEEN BASELINE AND OPTI- 
MIZED INLETS. 
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(g) X/D 2 = 3.3 


Mach: 0.02 0.15 0.29 0.42 0.56 0.69 


(h)X/D 2 =33 

External Surface Pressure 
Internal Surface Pressure 


FIGURE 13: MACH NUMBER CONTOURS AT VARIOUS 
CROSS SECTIONS, WITH SUPERIMPOSED CROSS FLOW 
STREAMLINES, IN EXTERNAL AND INTERNAL FLOW 
DOMAINS OF THE INLET. LEFT COLUMN: BASELINE IN- 
LET; RIGHT COLUMN: OPTIMIZED INLET. 


at about X/D 2 = 6.0, because the nacelle eventually completely 
submerges into the airframe and the external flow sees no ob- 
stacle, thereby returning to a flat surface and smooth flow. The 
exterior flow fields of both inlet are indistinguishable. 

On the other hand, the behavior inside the inlet is an en- 
tirely different matter. At the inlet entrance (X/D 2 = 0), the 
flow is also quite asymmetrical; but once inside the inlet, the 
flow appears to symmetrize relatively quickly within a short dis- 
tance, say roughly X/D 2 <0.5. However, a remnant of low-speed 
pocket from the entrance bottom wall first appears near the wind- 
ward side of the inlet, first inkling seen at X/D 2 = 0.5 in the 
optimized inlet and clearly visible at X/D 2 = 2.5 in both inlets; 
this low speed flow is of larger size in the optimized inlet. The 
growth of this momentum region persists, displaying spanwise 
nonuniformity; a secondary flow develops as the geometry un- 
dergoes a change from a noncircular shape to circular. Eventu- 
ally at AIP (X/D 2 = 3.3), the vortex becomes sufficiently strong 
that a vortex appears at the lower corner of the inlet, whereas 
a counter-rotating pair exist in the optimized inlet. Meanwhile, 
there appears a vortex adjacent to the sidewall on the leeward 
side. 

Additionally in Fig. [13] pressure distributions are plotted on 
the external and internal walls, denoted by black and orange col- 
ors respectively; it is noted that the high magnitude portion of 
the internal surface pressure turns out to belong to the upper sur- 
face. By and large both inlets behave similarly; some differences 
can still be discerned. At the entrance a noticeable difference 
in pressure is detected on the nacelle external surface, where the 
optimized inlet gives a relatively leveled distribution albeit some- 
what bumpy. At X/D 2 = 0.5, the lower surface of the optimized 
inlet also has more nonuniformity, explaining the earlier appear- 
ance of the low speed pocket mentioned above. 

Next we study the oil flow representations on floor surfaces 
in Fig. [l4]for both baseline and optimized inlets. Both are again 
similar, except inside the inlet entrance the optimized inlet ex- 
hibits spatial variations caused by modulating the wall shape. 
Both inlets reveal that the flow is significantly displaced outward 
from the nacelle in the windward side, resulting from the front 
flow separation discussed above regarding Fig. [13] The leeward 
side has got a large flow separation (X/D 2 ~ 2, see Fig. fl3l) , 
as a result the flow is greatly displaced away from the nacelle. 
There is a saddle-like singularity, located near the lip region re- 
spectively on the windward and leeward sides. The oil pattern 
on the bottom flow of the inlet is constricted considerably im- 
mediately after entering the inlet, showing a lateral push on the 
streamlines away from the wall. By comparing with the case of 
no cross wind, as given in Fig. \5\ we notice that the constriction at 
the low speed condition does not seem to be as severe inside the 
inlet as at high speed condition on cruise. Thus the total pressure 
recovery is higher at low speed. 

Finally, we show three views of oil flow patterns on the outer 
surface of the nacelle in Fig. [15] respectively: (1) the top view 
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p: 0.400 0.462 0.524 0.586 0.648 0.710 


FIGURE 14: OIL FLOW PATTERNS ON THE BOTTOM 
WALL OF INLET AND BOTTOM PLATE SUPERIMPOSED 
WITH PRESSURE CONTOURS. TOP: BASELINE INLET; 
BOTTOM: OPTIMIZED INLET. 


in the middle plot, (2) the sideview from the windward side in 
the top plot, and (3) the sideview from the leeward side in the 
bottom plot. A top view of the oil flow pattern on the plate is 
also included. The flow on the nacelle top surface behaves no- 
ticeably different from that on the inlet bottom wall (seen above 




p: 0.400 0.462 0.524 0.586 0.648 0.710 

FIGURE 15: OIL FLOW PATTERNS ON NACELLE EX- 
TERNAL SURFACES SUPERIMPOSED WITH PRESSURE 
CONTOURS. TOP: WINDWARD SIDE VIEW; MIDDLE: TOP 
VIEW; AND BOTTOM: LEEWARD SIDE VIEW. 


in Fig. IT4l) : the flow from windward side towards leeward side is 
seen, but the flow separation on the leeward side strongly pushes 
back, thereby showing converging of particle traces. Also, the 
two flow separation, each on the windward and leeward sides, 
shows the effects on the bending of particle traces. Similar be- 
havior is observed in the optimized inlet in these views, hence 
they are omitted. 
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CONCLUDING REMARKS 

Motivated by NASA’s interest in pursuing next generations 
of civil transport to environment and performance metrics, we 
have focused on the so-called N2B configuration which is a hy- 
brid (blended) wingbody with a top-mounted propulsion system 
into which a boundary layer of considerable thickness relative 
to the inlet height will be swallowed and diffused in an offset 
duct, resulting in a flow much more severely distorted than in the 
current fleet. This distortion results from a combination of fac- 
tors, such as vorticity embedded in the incoming boundary layer, 
lip separation, adverse pressure gradient, secondary flow, etc, 
all manifested into a enormously complicated flow phenomenon. 
Hence, it can be only analyzed via a high-fidelity setting to shed 
insight into the underlying flow physics. To keep flow distortion 
at the AIP within a tolerable level, some forms of control over the 
flow must be instituted, either by conventional (such as vortex 
generators, active plasma/synthetic jet, etc.) or other innovative 
ideas. Our previous studies show that it is incredibly effective to 
reduce flow distortion and increase total pressure recovery at the 
engine face by condi tioning/guiding the flow through modifica- 
tion of wall shape, instead of locally controlling through added 
penalizing devices. 

The objective of the present study is to assess the perfor- 
mance of the boundary-layer-ingesting inlet when subject to 
cross wind at low speed flight, e.g., during takeoff or landing. 
We focus on investigating the flow characteristics unique to this 
condition. Moreover, we are also interested in quantifying the 
characteristics in the inlet optimized for the cruise condition and 
comparing the differences in performance between the original 
baseline and the optimized inlets. In this investigation, we have 
shown that the flow in this off-design condition is very complex, 
involving additional flow separations external to nacelle, which 
will thus affect the aerodynamic performance of the aircraft to 
which the inlet is installed. As the flow enters the inlet with a 
crosswind component, its effect seems to last within a short dis- 
tance and quickly subsides by the constraining boundary of the 
inlet. At the low speed condition, the inlet distortion is low even 
with the cross wind effect, especially in the baseline inlet. The 
present study concludes that insofar as the inlet distortion is con- 
cerned, the crosswind effect, even up to 40 knots, is probably not 
of concern to operation of an engine fan. The pressure recovery 
in fact remains high, because of low losses owing to low speed 
flow. However, the distortion in the optimized inlet is several 
times higher than that of baseline inlet, although the pressure re- 
covery is only slightly lower. This suggests that a revisit to our 
previous optimization may be warranted; it is a follow-on inves- 
tigation that will be reported in a subsequent part of the series 
on the BLI inlet. In fact, this aligns with our view that for an 
effective mitigation of distortion in a BLI inlet system, a design 
not only considering the inlet passage itself, but also combining 
the airframe geometry surrounding the nacelle, is necessary. 
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